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O ' !• INTRODUCTION 

(N : 

Entanglement of quantum states is the most non-classical feature of quantum systems and one of the key resources in quantum 
information theory Q|]. In real quantum systems, inevitable interactions with surrounding environment may lead to decoherence 
\ resulting in degradation of entanglement. Since entanglement once it has been lost, cannot be restored by local operations, it 
is important to understand the process of disentanglement to control the effects of noise and to preserve as much entanglement 
as possible. The main motivation of the investigations presented in the paper is to study the possibility of the preservation of 
entanglement in the model of two two-level atoms interacting with environment with maximal noise (the case of thermal noise 
fNJ ' will be discussed elsewhere). In the Markovian approximation, the influence of that kind of environment on a single atom is 
described by dynamical semi-group with Lindblad generator in which the transition rates from ground state to excited state 
of the atom and vice versa, are equal to decay rate Fq. This noisy dynamics is related to the following limiting procedure: 
temperature of thermal reservoir tends to infinity, whereas the coupling strength goes to zero |2j |3[. On the other hand, the 
collective properties of two - atomic systems can alter the decay process compared with the single atom. It was already shown 
by Dicke [4] in the case of spontaneous emission and environment in the vacuum state that there are states with enhanced 
emission rates (superradiant states) and such that the emission rate is reduced (subradiant states). In the latter case, two-atom 
system can decohere slower compared with individual atoms and some amount of initial entanglement can be preserved or even 
, created by the indirect interaction between atoms 

In the model discussed in the paper, the collective dynamics of two atoms is described by damping rate F which depends 
on interatomic distance. When the atoms are separated by a large distance, one can assume that they are located inside two 
independent environments and T = 0. In that case, the noisy dynamics brings all initial states into unique asymptotic state which 
is maximally mixed. Moreover, all entangled states disentangle in finite time |9]. When atoms are confined in a region smaller 
than the radiation wave length, the collective damping rate is close to the decay rate, so we can use the approximation F = Fq. In 
that case, similarly as in the Dicke model, the antisymmetric state (singlet state) is decoupled from the environment and therefore 
is stable. This is the main physical reason why entanglement can be preserved in spite of the the influence of noisy environment. 

In this paper we are mainly interested in robust entanglement, so we study long time (asymptotic) behaviour of dynamical semi 
- group. We also not discuss here the case of arbitrary separation of atoms, since in that case all states disentangle asymptotically. 
Small distance separation modeled by the condition F = Fq, leads to the interesting semi - group which is not ergodic: asymptotic 
stationary states depend on initial conditions and can be parametrized by fidelity of initial state with respect to singlet state. We 
show that some of these asymptotic states are entangled. They belong to the important class of Werner states I1CB . In particular, 
we prove that if fidelity is greater then 1/2, the asymptotic Werner state is entangled. So collective damping can produce 
correlations between atoms which partially overcome the effect of decoherence, but in contrast to the zero temperature case 
HI IE Ulllll- this process cannot create entanglement. Our noisy dynamics has also another remarkable property: there are initial 
entangled states for which entanglement is preserved during the evolution, although the process of decoherence takes place, 
resulting in decreasing of purity. As we prove, asymptotic entanglement depends only on the overlap of the initial state with 
singlet state but not on its entanglement. So initial states with the same entanglement can behave differently with respect to the 
noise. This opens the possibility of protecting some entanglement of the initial state by performing local operations (which do 
not change entanglement) to maximize its overlap with singlet state. 

When the interatomic separation is small, coupling by the dipol - dipol interaction plays a significant role. It causes the 
entanglement between two atoms to oscillate in time, so for some period initial entanglement may even increase, but the noise 
decreases the amplitude of these oscillations, and asymptotically its contribution to the preservation of entanglement vanishes. 
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Since we mainly study here the role of noise in the time evolution of entanglement, we are not discussing in detail the interesting 
problem of the influence of dipol - dipol coupling on this evolution. Similarly, the experimental side of the problems studied 
here is beyond the scope of the present paper, where we investigate the theoretical model of the compound system interacting 
with maximally noised environment. 



II. NOISY DYNAMICS OF TWO-ATOMIC SYSTEM 



and 



where 



Time evolution of a density matrix of two two-level atoms A and B interacting with a noisy environment of the type discussed 
in the introduction, can be described by the following master equation 

^ = -i[H, p]+L N p (II.l) 

Here 

H = coo £ o ; 3 + £ Qjkoio- (H.2) 

j=A,B j,k=A,B,j^k 

l nP = ^ L ^jk(2aipa^-aiaip-paLa k + + 2a J _pa k + -a' + a k _p-pa J + a'L) (11.3) 

1 j,k=A,B 

o$. = o±®I, of=I<g>c ± , of = a 3 ®I, af = I<8>C3, a± = i (ai ± fo 2 ) 

and we identify ground state |0) and excited state |1) of the atom A or B with vectors ( ^ and f ^ in C 2 . In the hamiltonian 

term dll.21 . (Do is the frequency of the transition |0) — ► 1 1) and Q.ab — Qba = Q describes interatomic coupling by the dipol-dipol 
interaction. On the other hand, noisy dynamics is given by the generator dll.3t with 

Taa = r BB = r (n.4) 

and 

Tab = F ba =F (11.5) 

In the equality CI.4t . To is a decay rate of individual atom. The parameter F in the equality dll.51 describes the collective damping 
rate of two atoms interacting with a noisy environment. In our model, F satisfies 

r = g(R)F 

where g(R) is the function of the distance R between atoms such that g(R) — + 1 when R — ► 0. Notice that dll.3t can be rewritten 

as 

L N p- r [alpo* +a^po^ + a A p&l + a B pa B + -2p] 

+i r [2<po B + 2d pa B + - do^p - pdd - ddp - p<o B ] 

+ir [2dpd+2dpd-ddp-pdd-o^dp-po^d] (n.6) 

The master equation (III. It describing the time evolution of a two-atomic system in a noisy environment can be used to obtain 
the equations for matrix elements of any density matrix. To simplify calculations one can work in the basis of so called collective 
states in the Hilbert space C 4 1 1 1]. If 

/ 1 = |1)®|1),/ 2 = |1)®|0) ) / 3 = |0)®|1),/4 = |0}®|0) (II.7) 

then this basis containing excited state, ground state and symmetric and antisymmetric combination of the product states, is 
defined as follows 

\e) = /i » \g) = A, \s) = -j= {fz +h), \a) = ±= {h-h) (II.8) 
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In the basis of collective states, two-atom system can be treated as single four - level system with ground state \g), excited state 
\e) and two intermediate states \s) and \a) . From ( 111. It it follows that the matrix elements with respect to the basis \e) , \s) , |a) , \g) 
of the state p satisfy 

dpaa 



- (r - r) (2p aa - p gg - p ee ) 




(H.9a) 


-(T +r) {2p ss -pgg-p ee ) 




(H.9b) 


-iT p ee + (r + r) P , s + (r Q - r 


)Paa 


(H.9c) 


-2r pgg + (r + r)p ss + (r - r 


)Paa 


(H.9d) 


-(2r +4/coo)pe« 




(H.9e) 


-f2r n -2/£l)p„ s 

\ u /r US 




(II.9f) 


-(r -r)p, ?s + [(r+/i2)-2(r 


-i<Do)]pae 


(H.9g) 


-(r -r)p„ + [(r+in)-2(r 


+ i(O )]p ag 


(H.9h) 


(r +r)p,,-[(r+in)+2(r - 


2/COo) ] p re 


(H.9i) 


(r +r) P , 1 -[(r+/£2)+2(r + 


2/COq) ] p sg 


(n.9j) 



dt 

dpss 
dt 

dPee 

dt 

dpgg 

dt 

dp eg 

dt 

dPas 

dt 

dPae 

dt 

dPag 

dt 

dp se 

dt 
dPsg 

dt 

and for the remaining matrix elements one can use hermiticity of p. Notice that in jll.9t . equations flll.9al > - fll.9dt . ( I11.9el - 
dll.9fl and ( |II.9g| l - < |II . 9j | are decoupled and can be solved independently. Observe also that four - level system prepared in 
symmetric state \s) decays with enhanced rate To + T, whereas antisymmetric initial state \a) leads to reduced rate Tq — T. When 
two atoms are confined in a region smaller than the resonant wave length, we can put T = To (see e.g. (nj) so antisymmetric 
state | a) is completely decoupled from the environment and is decoherence - free state of the two - atomic system. 

One can also check that equations dll.91 describe two types of time evolution of the system, depending on the relations between 

r and To. When Y < To, there is a unique asymptotic state of the system, which is maximally mixed state In that case the 
relaxation process brings all initial states of two atoms into the state with maximal entropy. In the present note, we analyse the 
case when T = To and show that asymptotic stationary states depend on initial conditions and can be parametrized by matrix 
elements p aa of the initial states. 

m. SMALL DISTANCE SEPARATION AND GENERATION OF WERNER STATES 

When r = To, equations dll.91 simplify considerably and one can check that solutions of ( IH.9el > - ( |II.9j| > asymptotically 
vanish. Nontrivial contribution to the asymptotic stationary states comes from the matrix elements p„„, p„ , p ee and p gH satisfying 
equations dll.9ai - dll.9dl i. which in that case can be written as follows 



dp a 



dt 
dp^ 
dt 

dPee 

dt 
dpgg 
dt 



(III.l) 

= 2r (p <e + p w -2p M ) (111.2) 

2T Q ( Pss -pee) (1113) 

■ 2r (Pss-p gg ) (in.4) 
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The system of equations (IIII.H - dIII-4-i can be solved and we obtain 



Paa(t) 


= Pua(0) 




1 

— < 
3 


















p«(0 




P«a(0))H 


'- 6r °'(p«, 


(0) + 3p 4 


,(o)-i) 












Peeif) 




P«a(0))H 


1 


,-6r r (1 _ 


P«o(0) - 


3p*(0))H 


1 

~r 


- 2r °'(p flfl (o)- 


hp„(0)- 


l-2p ee (0) 


-1) 




-b- 


P«(0))H 


1 




p fla (0) - 


3p„(0))H 


i 

"2* 


" 2r °'(p««(o)- 


hp ss (0)- 


l-2p w (0) 


-1) 



(m.5) 
Cin.6) 

(m.7) 

(m.8) 



Thus 



limp.„(f) = -(l-p aa (0)) 

and similarly for p ee (t) and p g #(f). So the stationary asymptotic states p»o are parametrized by p au (0) = (a|p|a), where (a\p\a) = 
F is the fidelity of the initial state p with respect to the state \a) (or the overlap of p with singlet state \a)). In the canonical basis 
( III.7> . poo has the form 



l—F 














V o 



1 +2F 1-4F 

6 6 
1-4F 1+2F 



\ 
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1 —F 



(HI.9) 



Notice that for some values of parameter F, the state p„o is entangled. If we compute its concurrence given by the well known 
formula OEl 

C(p) = max (0, 2A, max (p) - trp ) 
where A, max (p ) is the maximal eigenvalue of p and 



P = V\/PPVP) P = (a 2 «>02)p(a2(8)a2) 
with p denoting complex conjugation of the matrix p, then we obtain 



C(Poo) 



0, F < 1/4 

2F-1, F>l/4 



So we see that for any initial state p with fidelity F, there exist asymptotic state p»o such that: 
1. if F £ [0, 1 /4), Poo is separable and can be written as 



(\ 











] P 2 P 



\ 






V o 



2 P 
T 



1-f 

o i+fy 



1 -4F 



(III. 10) 



2. ifF = |,p»=| 

3. if F £ (1/4, 1], poo is equal to the Werner state 



W a = (l-P)j+P\a){a\ 



(m. ii) 



with p = -^-^ — -. It is separable for F ^ (1/4, 1/2] and entangled for F > 1/2 with concurrence C(poo) = 2F — 1. 
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This result shows that all initially entangled states with fidelity greater then 1/2 preserve some amount of their entanglement 
during the interaction with maximally noisy environment and evolve into Werner states with the same fidelity. The class of 
Werner states has interesting properties: they interpolate between maximally entangled and maximally mixed states, for that 
class it was shown that some mixed entangled states can satisfy Bell inequalities 1 10], they have maximal possible entanglement 
with respect to the non - local unitary transformations and local and non - local general operations 1 14], Werner states can be 
also applied in entanglement teleportation via mixed states 1151 . 

Let us stress that the asymptotic behavior of the initial state depends only on its overlap with the singlet state \a) and not on its 
entanglement. There are many states with the same entanglement for which fidelity varies from to maximal value. This is for 
example the case of maximally entangled pure states. So initial states with the same entanglement can behave very differently 
with respect to the noise. The states which are more "similar" to the singlet state are more stable. Thus to protect the initial 
entanglement we may perform appropriate local operations on initial states to maximize its fidelity. All these aspects of noisy 
dynamics will be discussed on explicit examples in the next section. 



IV. SOME EXAMPLES 
A. Pure separable initial states 

Let p = |*F <g> <£>) (*F ® <E>| for l P, 4> € C 2 . Since for this state 

F = i(l-|(^,4>)| 2 ) (IV.l) 

< F < j and all asymptotic states are separable. Depending on the value of | (\P, <t>) | 2 we have the following possibilities: 

1. if i < |(\P, <E>)| 2 < 1, then the asymptotic state is equal to the state JIII.lQt with p = 2|( V P, <E>)| 2 - 1, 

2. if |(1',4>)| 2 = ^,thenpoc = ^ 

3. if < | (¥, 4>) | < then is equal to separable Werner state W a given by dill. 1 ll > with p = ^j- 3 — Notice that 

in contrast to the zero temperature case, where purely incoherent dissipative process can lead to the creation of entanglement 
El 01 [Hi in the present model initial separable states remain separable. 



B. Pure maximally entangled initial states 

Consider now the class of maximally entangled states flfill 

/ a 2 aVT^e-®' ay/T=2?e-*> - fl 2 e -'( e i+ e 2) \ 

1 



P(fl,9i,0 2 ) 



a a\ 
aVl-flV 9 ' 1 - a 2 

as 



-m 2 



iVY~a^e m ^ (l-a 2 )£>-'( e i- e 2) l-< 



V-flV'( e <+ e 2) -VT^e^ -aVT^e ie 



a 2 



(IV.2) 



where a € [0, 1], 9i , 02 € [0,2jt]. All states from the class (IIV.2I have concurrence equal to 1. On the other hand 

F = ^{\-a 2 ){\ -cos(6i -6 2 )) 



(IV.3) 



One can check that fidelity F can take all values from to 1 depending on parameters a and = 8i — 62- In particular F > j 
inside the set E on the (a, 9) - plane, given by 



1 a 2 a 2 

E = {(a, 9) : <a < —=, arccos < 9 < 2% — arccos -= } 

y/2 a — 1 a — 1 



Outside this set, F < j . On the curve 



2a 2 -1 



9 = arccos —-5 — — -, a e [0, V3/2] 
2\a L — 1 ) 



(IV.4) 



(IV.5) 



6>i — 02 




FIG. 1: Fidelity of maximally entangled states as function of a and 0i — 62. On the dotted curve F = ^ 

the fidelity is equal to | (see FIG.l). Thus all initial states P(a,0i,02) with (a, 61 — 62) <E E evolve into entangled Werner 
states, whereas states with (a, 81 — 82) outside E become separable. When these parameters lie on the curve dIV.5> . the dynamics 

brings corresponding initial maximally entangled states into maximally mixed state 5p For all initial states P(a, 81,82) with 
(a, 81 — 82) G E, asymptotic concurrence is smaller then 1, except antisymmetric state \a) which is stable. 

C. Some mixed initial states 

1. Mixed separable states 
If p is mixed separable state, then it can be written as 



P=L S *A' s k >0,Y,H = 1 

k k 



(IV.6) 



where are pure separable states. Since for any pure separable state, fidelity is not greater then 1/2, by jlV.6> the same is true 
for all mixed separable states. So they evolve into separable asymptotic states. 



Let 



Besides W a define also the states 



2. Werner states 



\±) = —( fl ±f 4 ) 



W s = (1-p)j + p \s) (s\, W± = (l-p)j+p \±) (±| 



(IV.7) 



One can check that for all states flIV.7l > 
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so < F < \ and they evolve to asymptotic separable state JIII. 10i . On the other hand, the states (IIV.7I are locally equivalent to 
W a . If we define 

U s = 03(g>I 2 , U + = I 2 ®ics2, C/_=l2®ai 

then 

W s = U s W a U*, W+ = U+W a U*, W- = U-W a U*_ 

3. Bell-diagonal states 

Let ps be the convex combination of pure states |+), |— ), \s) and \a) 

Pb = Pl\+)(+\ +P2|-X-| +P3\s)(s\ +p 4 \a)(a\ (IV.8) 

It is known that if all p,- € [0, 1/2], Pn is separable, while for p\ > 1/2, is entangled with concurrence equal to 2p\ — 1 
(similarly for P2, P3, Pa) 11*711 . On the other hand, for states JIV.8I 

F = P4 

so all states (separable or entangled) with p4 < 1 /2 become separable asymptotically. If p\ > 1/2 then noisy dynamics produces 
asymptotic Werner state W a with concurrence equal to 2F — 1 = 2p4 — 1. In this case asymptotic entanglement is exactly equal 
to initial entanglement, so the amount of entanglement is preserved. In the next section, we discuss this interesting phenomenon 
for some class of initial states. 

V. PRESERVATION OF ENTANGLEMENT FOR SOME NON-MAXIMALLY ENTANGLED INITIAL STATES 

Let us discuss now how entanglement of the asymptotic state can depend on initial entanglement. This problem has a simple 
solution in the case when initial entanglement is a function of fidelity F. Since F is constant during the evolution, the final 
entanglement is exactly equal to its initial value. So the process of collective damping can preserve entanglement of some initial 
states. Simple examples of such states are described below. Let 



/0 0\ 

P22 P23 

P23 P33 

\0 0/ 



(V.l) 



Notice that in JV. li all matrix elements are real and F = 1(1 — 2p23). Moreover 

C(p)=2|p23| = |l-2F| (V.2) 

From our previous results it follows that: 

a. if p 2 3 > i.e. F < \, then C(p) = 2p 23 > 0, but C(poo) = 0, 

b. if p 2 3 < i.e. F > \, then C(p) = -2p 23 = 2F - 1 and C(p„) = 2F - 1, so 

C(p.)=C(p) (V.3) 



Notice that the unitary operator Oj, ® I2 transforms the states JV. li with P23 > to the states with P23 < 0. So performing local 
operations on the initial state, we can protect its entanglement. 

Consider also explicit examples of states for which the relation dV.31 holds. Let us take the class of pure states 

l P = cos^|0)(8)|l>+sin^|l)(8)|0), 0G[O,Jt] (V.4) 

One can check that for the class dV.4l > 

F = i(l-sin24») 
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and 



C(poc) = max (0, — sin2(])) 

We see that all pure entangled states (IV. 4> with <|) 6 [7t/2,7t], evolve into asymptotic Werner states, which have the same entan- 
glement as initial states. 

It is instructive to discuss evolution of initial states iY.lt also for finite times. Here the dipol - dipol interaction between atoms 
plays a significant role. For small separation of atoms the ratio £2/To can be large 1 1 ll fl8fl and this coupling can induce transient 
increase (or decrease) of entanglement. But noise reduces the amplitude of those oscillations of concurrence and they vanish 
asymptotically. To see some details, we solve the equations IIII. li — ( 1111. 4> and ill. 9ft for such initial states and obtain 



p(0 



/Pn(f) o \ 

P22(f) P23(f) 

P32(f) P33(0 

V o o o P44 (f)/ 



(V.5) 



where 



pn(0 


= i(l+2p 23 ) 


1 

~v 


; - 6r «'(i 


+ 2p 23 ) 






P22(?) 


= i(l-p 23 )H 


i 

— e 
6 


- 6r »'(H 


-2p 23 )- 


f -e" 2Vot cos2£2f(p 22 - 


P33) 


P33(0 


= |(1-P23)H 


1 


-6Tof( l4 


■2p 23 )- 


-i e - 2r ° r cos2£2f(p 22 - 


P33) 


p44(0 


= i(l+2p 23 ) 


i 

-? 


r 6r 0<(1 


+ 2p 23 ) 






P23(0 


= i(4p23"l) 


i 

+ r 


>- 6r °'(l 


+ 2p 23 ) 


_i e -2rv sin2 £2f(p 22 - 


-P33 



One can check that for states dV.5> . concurrence is given by 

C(p(0) - max (0, 2 ( |p 23 (f)| - VPn(0p44(0 



so 



with 



= max(0,2(|p 23 (f)| -pii(f))) 
C(p(f))=max ^0, 2y / A(f) 2 +B(f)2 s in 2 2af-2p 1 i(f)^ 



A (t) = l (4p 2 3-l+e- 6r< "(l+2p 23 ) 



B{t) = -e- 2T « (p 22 -p33) 



(V.6a) 
(V.6b) 
(V.6c) 
(V.6d) 
(V.6e) 

(V.7) 
(V.8) 

(V.9) 



This shows that the concurrence as a function of time oscillates with period depending on the strength of dipol - dipol interaction 
but when p 23 < 0, it tends to its initial value (FIG. 2). When the dipol - dipol interaction is absent or when the initial state is 
such that p 22 = P33 = j and p 23 < 0, then from lV.8t and i\.9l it follows that 



Ip23(0|-Pn(0 = -P23(0-Pn(0 



(V.10) 



Since time-dependent terms in ( IV. lOi cancel, it is equal to — 2p 23 , and 

c(p(0) = c(p) 

for all t . When p 23 > 0, then 

|P23(0|-PH(0=P23(0-PH(0 



(V.ll) 



and ( IV. lll idepends on t in such a way that C(p(f)) monotonically goes to zero. The effect of preservation of entanglement for 
states dV.5l > should be contrasted with the monotonic decreasing of purity of the state p defined by trp 2 , during the time evolution 
( III. II . As shown in 1 19], it is equivalent to the condition 

L N (I 4 ) = 

which can be simply checked to be true in our model. 
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0.25 0.5 0.75 1 1.25 1.5 1.75 2 1 0^ 

FIG. 2: Concurrence as function of time for initial state <V.4> with = 27t/3 for fl/To = 1 (solid line) and fl/To = 3 (dotted line) 

VI. CONCLUSIONS 

During the evolution of open quantum systems interacting with environment the process of decoherence usually results in 
degradation of entanglement. To preserve as much entanglement as possible one has to control the effects of noise. In this 
context, we have studied the model of compound system of two atoms influenced by so called maximal noise, which can be 
treated as the limiting case of thermal noise, when temperature goes to infinity. As we have shown, even in that case there are 
entangled states which are decoherence - free. Explicit examples are given by singlet state \a) and Werner state W a . On the other 
hand, there are evolving states with a very interesting property: its asymptotic entanglement is exactly equal to the initial one, 
or even is stable during the time evolution. We have also shown that performing some local operations on initial states can help 
with protecting entanglement. 
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